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ABSTRACT 



Context. A bright feature 80 pc away from the core in the powerful jet of M 87 shows highly unusual properties. Earlier radio, optical 
and X-ray observations have shown that this feature, labeled HST-1, is superluminal, and is possibly connected with the TeV flare 
detected by HESS in 2005. It has been claimed that this feature might have a blazar nature, due to these properties. 
Aims. To examine the possible blazar-like nature of HST-1, we analyzed A2 cm VLBA archival data from dedicated full-track obser- 
vations and the 2 cm survey/MOJAVE VLBI monitoring programs obtained between 2000 and 2009. 
Methods. Applying VLBI wide-field imaging techniques, the HST-1 region was imaged at milliarcsecond resolution. 
Results. Here we present the first 2 cm VLBI detection of this feature in observations from early 2003 to early 2007, and analyze its 
evolution over this time. Using the detections of HST-1, we find that the projected apparent speed is 0.61 ±0.3 Ic. A comparison of the 
VLA and VLBA flux densities of this feature indicate that is mostly resolved on molliarcsecond scales. This feature is optically thin 

(a 0.8 for S oc v*") between A2 cm and A20 cm. 

Conclusions. We do not find evidence of a blazar nature for HST-1. 
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1. Introduction 

Active Galactic Nuclei (AGN) are among the most ener- 
getic phenomena in the Universe, and they have been heav - 
ily studied since their initial discovery by ISevfertl (Il943h . 
Although there are many clues that imply that a super-massive 
black hole (SMBH) is t he engine launching the powerful jet 
(lUrrv & Padovanil Il995h . the exact mechanism remains un- 
known. 

One of the most studied AGN is M87 (also known as 
Virgo A), a nearby elliptical galaxy located in the Virgo clus- 
ter. It hosts a very powerful one-sided jet emerging from the 
central r egion, and wa s also the first extragalactic jet to be dis- 
covered (ICurtisl[T918h . T he synchrotro n-emitting nature of the 
M 87 jet was suggested bv lBaadel(ll956h . Observations show that 
M87 contains 2.4 x 10^ Mq within a 0.25" radius, which sug- 
gests the presence of a SMBH in its nucleus dHarms et al. 1994). 
Due to its proxim ity (16.4Mpc, z= 0.00436, 1 mas = 0.08 pc, 
1 mas yr -1 = 0.26 c, [Jordan et al.ll2005 ). M 87 is an ideal candi- 
date for studying AGN phenomena, and has been monitored at 
many different wavelengths over the last several decades. In the 
observed one-sided jet of M87 (Shklovskv 1964), superlumi- 
nal motion was reported from Hubble Space Telescope (HST) 
observations within the innermost 6" of the jet with apparent 
speeds of 4 c to 6 c dBiretta et alj|1999l) . Discrepant speeds were 
reported from VLBA observa tions at 7 mm with values between 
0.25 c to 0.4 c dLv et alJ l20()7h. and a value of 2 c or even larger 
dAcciari et alj|2009h . VLBA /}2cm observ ations showed appar- 
ent speeds < 0.05 c from 1994 to 2007 dKovalev et alJ l2007t 



iLister et al]l2009bl) . Therefore, the kinematical properties of the 
jet in M 87 remain an important topic of discussion. 

In 1999 HST observations revealed a bright knot in the jet 
located 1" (projected distance of 0.08 kpc) away from the core. 
This feature, named HST-1, is active in the radio, optical and X- 
ray regimes. VLBA A2Q cm observations showed that HST-1 has 
sub-structure and appears t o contain sup erluminal components 
moving at speeds up to 4 c (Cheung et alj|2b07i) . These observa- 
tions suggest that HST-1 is a collimated shock in the AGN jet. 
Furthermore, recent multi-wavelength observations were used 
to show that HST-1 could be related to the origin of the TeV 
emission observed in M 87 in 2005 by the HESS telescope 
( Aharonian et al. 2006). Comparing data taken in the near ultra- 
violet by Madrid (2009), soft X-rays (Chandra) imd V LA A2 cm 
observations dCheung et al.ll2007t [Harris et al.ir2006h . the light 
curves of HST- 1 reached a maximum in 2005, while the resolved 
core showed no correlation with the TeV flare. Therefore, the 
TeV emission from M 87 was suggested to originate in HST-1 
dHarris et al.ll2008h . Based on those findings. lHarris et al.1 d2008h 
sugges t that HST-1 has a blazar nature. However. lAcciari et al.l 
(2009') reported rapid TeV flares from M 87 in February 2008, 
which may have originated in the core, instead of HST-1, which 
remained in a low state during the flares in 2008. Recently, the 
Fermi Large Area Telescope (L AT) team reported the detection 
of M 87 at gamma-ray energies dAbdo et alJl2009l) . 

The AGN standard model considers the blazar behavior to 
originate at the vicinity of the SMBH. However, HST-1 is 80 pc 
away from the core. If the HST-1 blazar hypothesis is true, this 
would pose a challenge to current AGN models. In this paper. 
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Table 1. Journal of VLBA A2 cm observations of M 87'' . 
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" Observations are in dual-polarization mode, except as indicated 

Total scheduled VLBI on-source time 
^ Beam A: natural weighting and tapering of Gaussian full-width half maximum 0.3 at 200 M/i 

Root-mean-square image noise using natural weighting and tapering with a Gaussian factor of 0.3 at a radius of 200 M/l 

Recording rate: 256 Mbit s"', 2 bits per sample 

' Dedicated full-track experiment on M 87 (see lKovalev et al.ll2007h : the array included a single VLA antenna (Yl) 
^ Recording rate: 128 Mbit s"' , 1 bit per sample 

Single polarization hand (LL) observations only 
' Recording rate: 512 Mbits"', 2 bits per sample 
j Antenna NL missing 
^ Antenna SC missing 
' Antenna HN missing 



we examine this hypothesis with A2 cm VLBI wide-field imag- 
ing of the HST-1 feature. We present the description of our VLBI 
observations and the corresponding data reduction in Sect. |2j 
the results are presented in Sect. |3j the discussion is presented 
in Sect. |4] Finally, a short summary is presented in Sect. |5] 
Throughout this paper, we use the term "core" as the apparent 
origin of AGN jets that commonly ap pears as the brightest fea- 
ture i n VLBI radio images of blazars jLobanovIl 1 998tlMarscherl 
120081) . We use a cosm ology with Q„, ^ 0.2 7, Qa = 0.73, and 
Ho = 71kms-i Mpc-i (iKomatsu et al.ll2009l) . 

2. VLBI Observations and Data Analysis 

M 87 has been monitored at A2 cm with the VLBA since 1994 by 
the 2 cm Survey /MOJAVE program^] (Kellermann et al. 20o4 
iLister et al. 2009a). We re-analyzed 12 epochs of these monitor- 
ing program data obtained after late 2001, together with three 
observing sets of targeted observations on M87 in 2000 (see 
Table©. The 2 cm Survey/MOJAVE epochs from 2001 to 2009 
were snapshot observations with a total integration time from 20 
minutes to one hour, and were broken up into about six-minute- 
long scans in order to have a better (m, v) coverage; the three 
epochs in 2000 were full-track observations with eight-hour in- 
tegration time and included a single VLA antenna (Yl). 

The VLBA data were processed by following the stan- 
dard procedures in AIPS cookbook, as described in detail by 
iKovalev et al](l2007l) : lLister et al.l ( l2009ah . The data were fringe- 
fitted before the imaging process. The data from MOJAVE 
project epochs after 2007.61 were processed using the pulse 
cahbration signals to align the phases instead of fringe fitting, 

' |http : //www ■ physics . purdue . edu/MOJ AVE/ 1 



because the position s of the sources and the VLBA antennas 
are well determined dPetrov et al.ll2008l l2009l) , and the station 
clocks are known to be stable. However, for our wide-field imag- 
ing purposes, fringe fitting was necessary for a better determina- 
tion of the phase and rate across the observed band, and therefore 
to better image the extended jet of M 87. Moreover, it ensures the 
homogeneity of all the datasets. 

In order to accurately image HST-1 with VLBI, there are two 
issues that need to be considered. First, HST-1 lies ~800 A2cm 
beamwidth away from the brightest feature (the VLBA core). 
For this reason, time or frequency data averaging would pro- 
duce time and bandwidth smearing in the HST-1 region. Second, 
based o n earlier VLBA A2Q cm and VLA A2 cm observations of 
HST-1 (ICheung et al.ll2007h . we expect the total 2 cm flux den- 
sity to be at the milliJansky level, which is much weaker than 
the total flux density of the inner jet (~2.5 Jy). To detect HST-1, 
we need to image the inner jet region with its extended structure 
(over tens of milliarcseconds); otherwise, the sidelobes from the 
core would cover the HST-1 emission. To reach this goal, we ap- 
plied natural weighting and (m, v)-tapering on the whole dataset. 

First, we averaged the data 16-sec in time and averaged all 
channels in each IF in frequency, and obtained the inner jet ini- 
tial model by applying phase and amplitude self-calibration us- 
ing DIFMAP. The time-smearing effect of 16 seconds begins 
to be relevant at a distance of 200 mas from the field center, 
which is beyond the inner jet scale, and did not affect our data 
analysis at this stage. When the inner jet model was reasonably 
good, we applied the obtained CLEAN-components to the un- 
averaged data for the first amplitude-and-phase self-calibration 
using AIPS. Natural weighting and tapering with a Gaussian fac- 
tor of 0.3 at a radius of 200 MA in the (u, v)-plane was applied 
to the whole dataset, and the resultant beam sizes are shown in 
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Fig. 1. Images of the central 220 mas of the M87 jet (beam A, 
resulting from a [u, v)-taper of Gaussian 0.3 at 200 M/l and nat- 
ural weighting), spaced by their relative time intervals. The im- 
ages are plotted with the same size scale. The contour levels all 
increase by factors of V2, and the lowest contours are (from 
2000.06 to 2008.33): 0.5, 0.7, 0.5, 1.1, 1.3, 0.9, 0.9, 0.9, 0.9, 
0.9, 1.1, 1.0, 1.2, 0.9, 0.9 mJy beam"', respectively. The restor- 
ing beams and brightness peak values are given in Tabl^U 



Fig. 2. Downgraded resolution images of the central 220 mas 
of the M87 jet with beam B (8x3.4 mas), spaced by their rel- 
ative time intervals. The contour levels all increase by factors 
of V2, and the lowest contours are (from 2000.06 to 2008.33): 
0.5, 0.7, 0.5, 1.1, 1.3, 0.9, 0.9, 0.9, 0.9, 0.9, 1.1, 1.0, 1.2, 0.9, 
0.9 mJy beam"', respectively. 



Table [T] If we remove tapering and apply uniform weighting, the 
inner jet is better resolved but less flux density is recovered. In 
this case, the image noise level is too high to detect the HST-1 
feature. 



To better identify the extended component, and since full- 
resolution, uniform-weighting imaging did not produce detec- 
tions of this feature, we chose the approach of tapering and nat- 
ural weighting. Furthermore, we downgraded the resolution to 
a larger beam size of 8x3.4 mas in P. A. 0° This beam size was 
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Fig. 3. Total VLBA flux density of M87 at XI cm versus time. 



chosen specificall y for conipariso n with the VLBA /I20cm ob- 
servations of (Cheung et alJl2007l) In the following Sections, we 
label the smaller beams as beam A, noting that the beam sizes 
are slightly different in each epoch. These are listed in Table [1] 
Beam B refers to the beam size of 8x3.4 mas. 



3. Results 

We completed the wide-field imaging of 15 epochs of M87 
VLBA XI cm obtained between 2000 and 2009. Each epoch 
was imaged using beam A and beam B (see Sect. |2]i, and for 
each beam size, we applied two image cleaning fields, one for 
the M87 inner jet region, and another for the HST-1 region, 
which was phase-shifted as -788.5 mas in right ascension and 
348.9 mas in declination. We performed deep cleaning iterations 
until the image rms reached the expected thermal noise value. 
The rms values of the final images are listed in Table [1] Three 
epochs in the year 2000 are full-track observations with an on- 
source time of 476 minutes, and have the lowest rms. 



3. ^ . Imaging of the inner jet of M87 

With tapering and natural weighting of the data, the extended 
inner jet was imaged over a region 100-200 mas in size, which 
varied according to the different sensitivity of each epoch. The 
inner jet images are shown in Fig. [1] (beam A) and Fig. |2] (beam 
B). The inner jet structure has an average total flux density of 
~2.5 Jy, with overall flux density changes up to 1 Jy during 2000 
to 2009, as shown in Fig.[3] We estimate the error in flux densi- 
ties to be 5%, base d on the typical am plitude calibratio n accu- 
racy of the VLBA ( Kovale v et ai]|2005l and|Ulyestad & Wrobell 
12009 ). In February 2008, multiple very high energy flares were 
detected by the multi-wavelength campaign of HESS, MAGIC, 
VERITAS; meanwhile, VLBA XI mm observa tions detected a 
flare in the VLBA core simultaneously dAcciar i et al. 2009). We 
have also seen the 2008 flare from the core in our VLBA A2 cm 
data, and therefore confirmed this result (see Fig.[3]|. 

The inner jet showed structural changes during the period 
of our observations. However, the study of the inner jet is be- 
yond the sc ope of thi s paper (see e.g., Kovalev et al. 2007. 
I Walker et aLl l2008, and lLister et al. 2009b for a detailed discus- 
sion of the structure and kinematics of the inner jet). 



3.2. Imaging of the HST-I region 

HST-1 was detected in six of the 15 epochs we analyzed from 
2003 to early 2007 (detection limit: 5cr). As discussed in Sect.|4] 
this feature was too faint in most epochs, and the image noise 
was too high to allow us to to detect it. Fig.|4] shows hybrid 
maps of this region. The HST-1 images produced using beam A 
show that the brightest component has a size of ~ 10 mas down 
to the lowest contours, while the images produced using beam 
B reveal an extended structure ~50mas in size, which is com - 
parable with VLBA /120cm observations (Che ung et al.ll2007h . 
The feature has a total flux density that varies between 4 and 
24mJy (Fig.|5]l. The peak surface brightness varies from 1 to 
4 mJybeam"' (beam A) and 2 to 10 mJybeam"' (beam B), as 
shown in Fig.|6] Epochs which have no detection are marked as 
upper limits (inverted triangles) in Figs.|5]and|6] 

In 2003.09, HST-1 is marginally detected with both beam 
sizes; the brightness peak is ~1 mJybeam"', and the total flux 
density is 3.7 mJy. During 2004.61-2005.85, 4 epochs show sig- 
nificant detections of HST-1 with beam A and beam B, where 
the peak component remains dominant, and the total flux den- 
sity fluctuation is -30%. In 2007, the structure of HST-1 be- 
comes more extended and complex, reaching a total flux density 
of 10 mJy and developing two distinct peaks (peak a and peak b 
in Fig. nil. 

For the epochs with no detections, we derived upper limit 
on the total VLBI flux density of HST-1. These were estimated 
based on the rms noise values of each epoch. The maximum 
amount of flux density hidden below the noise would be the rms 
value times the size of the HST-1 emission region (in units of 
beam size). Based on our HST-1 detection, we assume that the 
size of the HST-1 emission region is about 20 times the beam 
size, and we derived the upper limit on each epoch accordingly. 
We also derived the upper limit of the brightness temperature in 
this region as 9x10'' K at X2cm. 



3.3. Spectral properties of HST-I 

Figurelll illustrates two overlaid images of HST-1 at XZQcm 
(Cheung et al. 2007 and Cheung, priv. comm.) and A2 cm (beam 
A) of two adjacent epochs. The images were co-aligned using 
the peak of the inner jet as a common reference point. This figure 
shows that our observations are resolving the HST-1 peak com- 
ponent. We use those epochs to derive the HST-1 spectral proper- 
ties by producing a rendition of the A2 cm image using the same 
restoring beam as the A20 cm image, namely, 8x3.4 mas at a po- 
sition angle of 0° (beam B; Fig.|4] right panel). Table|2]lists the 
epochs of A2Q cm and A2 cm observations and the coiTespond- 
ing spectral index a, where SvOcv+". The uncertainties of the 
spectral indices are formal errors, and are estimated using stan- 
dard error propagation methods. When interpreting these values, 
one has to keep in mind that the resolution of A20 cm and A2 cm 
are very different. First, the frequency is different by one mag- 
nitude, which results in small formal eiTor bars of the spectral 
index. Second, the incomplete (m, v) coverage on short VLBA 
baselines might cause a missing flux density in the HST-1 re- 
gion at A2 cm (see Sect. l4.4l for a discussion). The average value 
of the derived spectral index is a 0.78. 



3.4. HST-I kinematics 

To study the kinematics of HST-1, we fitted the peak of HST-1 
and the M 87 core in the image plane with a Gaussian component 
using IMFIT in AIPS. As shown in Fig.H] it is difficult to iden- 
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Fig. 4. Images of the HST- 1 region restored with beam A (left panel) and beam B (right panel). The distance between epochs 
is proportional to the relative time interval, and the images are plotted with the same size scale. The contour levels increase by 
successive factors of V2, and the lowest contours for both images are (2003.09 to 2007.10): 0.6, 0.8, 0.8, 0.8, 0.8, 0.7 mJybeam"\ 
respectively. The beam sizes of the beam A images are plotted to the right of each epoch label, while the beam size of the beam 
B images is shown in the upper-right corner. The HST-1 field is phase-shifted -788.5 mas in right ascension (RA) and 348.9 mas 
in declination (DEC). These images were obtained after performing self-calibration in the imaging process. In the right panel, a 
and b represent the two peaks of HST-1 in epoch 2007.10. The red dashed lines illustrates the linear fit yielding HST-l's projected 
apparent speed (see Fig.|8]i. 
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Table 2. HST-1 ^20-/12 cm spectral index^ 



Epoch 


Sv [mjy] 


a 


A2 cm A20 cm 


A2 cm A20 cm 




2005.30 2005.35 


19.7±1.0 111±6'' 


-0.75±0.03 


2005.85 2005.82 


19.9±1.0 126±6'' 


-0.80±0.02 



The spectral index and its errors are formal, see section 14141 
5% error assumed for VLBI observations 
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Fig. 5. Total VLBI A2 cm flux density of the HST-1 region versus 
time. The values for epochs with no confident detection (inverted 
triangles) are upper limits, whereas the epochs with a detection 
are marked as dots. The results of using beam A are in black 
closed circles/triangles, while the results of using beam B are in 
red open circles/triangles. 



tify moving components in this region, because the detections 
are weak. Therefore it is difficult to derive an accurate apparent 
speed for the HST-1 subcomponents. For this reason, we use the 
relative position between the peak of HST-1 and the M 87 core 
to estimate the apparent speed of HST-1. Moreover, among the 6 
epochs of detection, HST-1 was almost resolved out in epoch 
2007.10, and appeared to have a double-peak morphology in 
beam B image (Fig.|4]i. By considering that the HST-1 detec- 
tion in 2007.10 was weak, and the double-peak structure might 
be due to the effect of convolution, we excluded this epoch from 
the kinematic analysis (see Sect. l4.2l for a discussion). As shown 
in Fig. [8] we plotted the peak positions of HST- 1 with respect to 
the M 87 core against time, and we estimated the position errors 
as the FWHM of the fitted Gaussian component. By applying a 
Unear regression to the peak positions, we obtained a value for 
the projected apparent speed of the HST-1 feature of 2.36±1.19 
masyr"', which correspond tOjSapp = 0.61+0.31. 
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Fig. 6. Brightness peak value at VLBI A2 cm of the HST-1 region 
versus time. The values for epochs with no confident detection 
(inverted triangles) are upper limits, whereas the epochs with a 
detection are marked as dots. The results of using beam A are in 
black closed circles/triangles, while the results of using beam B 
are in red open circles/triangles. 



size where the flare was produced. We define the logarithmic 
variabiUty timescale as Tvar - dt/d[ln{S)], where S is the flux 
density, t is the time interval between observations in units of 
years (Burbidgeet al. 1974). Next, we estimate the upper limit 
of the characteristic size of the emission region from light-travel 
time: flchar ~ 0.13 (1 + z)D 'Tvar<J, where z is the redshift, D 
is the luminosity dista nce in Gpc, and 6 is the Doppler factor 
dMarscher et al.l ll979l) . Table[3] shows the result of the derived 
ch aracteristic si z e of th e 2005 flare in HST-1. 

[Harris et all ( l2003h estimated the Doppler factor 6 of HST- 
1 to have a value between 2 and 5, based on de cay timescales 
of Chandra observations. In the same context, I Wang & Zhoul 
(2009) estimated the Doppler factor of HST-1 to be 3.57±0.51 by 
fitting the non-simultaneous spectral energy distribution of M87 
using a synchrotron spectrum model. Simulations incorporating 
MHD models for the M 87 jet have suggested that the recollima- 
tion shock formed close to the H ST-1 position had a relatively 
low Doppler factor -1-2 dOracia et al.ll2009 ). If we use (5=3.57, 
the derived characteristic sizes of HST-1 emission region during 
its flaring time are 20 < On^e <31 and 54 < < 100 mas. The 
size scale of structural changes of HST-1 (see Fig.|4|l is between 
20-50 mas (1 pc = 12.5 mas), which is within the derived char- 
acteristic size scale. However, if 6 <3.5, ^rise derived from VLBA 
/12 cm is smaller than 20 mas, which would create a causality 
problem, since the largest structural change of HST-1 cannot be 
bigger than the upper limit of the information propagation time. 
Therefore, 6 >3.5 is needed based on the causality argument. 



4. Discussion 

4.1. Variability timescale and HST-1 flaring region at 
parsec-scales 

We derived the variability timescale of the 2005 flare from HST- 
1 assuming that a single flare would produce a logarithmic rise 
and fall in the light curve. By fitting the available light curves 
from VLBI and VLA measurements before and after the max- 
imum flux density during the flare, we derived corresponding 
variability timescales and the upper limit of the characteristic 



Table 3. HST-1 variability timescale and characteristic size dur- 
ing the rise and fall of the 2005 flare. 



Wavelength 




Tfall 








[yr] 


[yi] 


[mas] 


[mas] 


VLA A2 cm 


1.5 


4.9 


8.7(5" 


28(5 


VLBA A2cm (beam A) 


1.0 


4.2 


5.6(5 


24(5 


VLBA A2cm (beam B) 


1.0 


2.7 


5.7(5 


15(5 



"The characteristic size is related to the Doppler factor 6. 
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Fig. 7. VLBA images of the HST-1 region in M 87 at /1 2 cm (this paper: black contour, beam size at half power level = 2x1 mas, 
P. A. = -16°, plotted bottom left) and /120cm (ICheung et alJ (2.00 7): gray contour, beam size = 8x3.4 mas, P. A. =0°, plotted bot- 
tom right). Left-hand panel: epoch 2005.30 (/I2cm, peak surface brightness: 3.4mJybeam ') and 2005.35 (/i20cm, peak: 
45 mJy beam"'); right-hand panel: epoch 2005.85 (/12cm, peak: 3.3 mJy beam"') and 2005.82 (/i20cm, peak: 42mJybeam"'). 
The lowest contour is 0.7mJybeam"', and the contour levels are separated by a factor of V2. The images of two frequencies were 
registered using the peak of the inner jet of HST-1 as a common reference point. 



If 6 is less than 3.5 in the HST-1 region, we could explain 
the causality problem by the self-calibration procedure that we 
used while imaging HST-1. We applied self-calibration to the in- 
ner jet and HST-1 for all of the epochs with detections in beam 
A and beam B. Self-calibration works better with stronger ob- 
jects, in our case, the marginal detection in epoch 2003.09 might 
not provide a satisfactory result in recovering extended emis- 
sions, comparing with stronger detections between 2004.61 and 
2005.85, in which after self-calibration, we were able to recover 
more extended emission. 



Id has jSapp = 1.14 + 0.14. Our results are consistent with these 
findings. 

4.3. Detection limits 

HST-1 was not detected in epoch 2006.45. This epoch had 
only 22 minute integration time and the sampling rate was 
128Mbit s"' . Therefore, the rms level (0.35 mJy beam"') was not 



4.2. Speeds of HST-1 

We used the peak position and fitted Gaussian errors to estimate 
the projected apparent speed of HST-1 jS-^pp = 0.61+0.31 (see Fig. 
[8]l, which is sub-luminal. In the kinematic analysis, we excluded 
epoch 2007.10 for obtaining a robust kinematic result. However, 
to make our results more solid, we estimated the upper and lower 
limit of HST-1 apparent speed by including epoch 2007.10 as a 
test. The derived possible range of the apparent speed is 0.23 



< A 

jet motion. 



app 



1 .2, which is still consistent with a mildly relativistic 



To compare our kinematic results with previous findings, 
we show the positional evolution on the sky of the HST- 
1 peak at VLBA A2cm with both beam A and beam B 
(see Fig. [Hi, and c omponents CI and C2 at VLBA /I20cm 
dcheungetalJ 120071 and priv. comm.). As illustrated, our de- 
rived apparent speed range and structural evolution in time at 
A2 cm are consistent with the /120 cm results. However, our sub- 
luminal speed measurements of HST-1 are inconsi stent with the 
high apparent superluminal motions reported in ICheung et al.l 
(2007) from VLBA /120 cm observations and with the HST 
(Biretta et al. 199 9). Nevertheless, in the VLBA /120cm obser- 
vations dCheung e t al. 2007), lower speeds were derived from 
some components in HST-1: c2 has jSapp = 0.47 ± 0.39, and HST- 



_ HST— 1 Kinematics 



/^app = 0.61±0.31 



CS870 



O 

s 



850 



2003 2004 2005 2006 2007 

Time (yr) 

Fig. 8. The linear fit of HST-1 proper motion shown as the red 
dashed line. This plot illustrates the position of HST-1 compo- 
nent peaks with respect to the M 87 core component as a function 
of time from 2003.09 to 2005.85. The fitted projected apparent 
speed of HST-1 jSapp is shown at the upper- left corner of the plot. 
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Fig. 9. HST-1 sky p ositions (left) and radial distance as a function of time (right) relative to the M87 core. VLBA /I20cm 
dCheung etal J 120071 priv. comm.): time range of components CI and C2 is from 2005.04 to 2006.53; VLBA /I2cm (beam B): 
time range of components a and b is from 2003.09 to 2007.10; components a and b are shown in Fig.H] 



low enough to detect HST-1, which was also fading according to 
our light curve (Fig.|5]). 

The total flux density of HST-1 measured by the VLA 
at /t2cm rea ched its maximum value of ~123mJy in 2005 
dHarris et alJ l2006). Our VLBA A2 cm results recovered a HST- 
1 total flux density of ~23 mJy, which is only 19% of the VLA 
measurement. We conclude that the innermost region of HST-1 
was resolved out with the long baselines of the VLBA. The low 
measured flux density suggests that HST-1 is very extended. 

4.4. HST-1 parsec-scale spectrum 

We derived the spectral index of the HST-1 region based on 
VLBA A2cm and /I20cm data (Table |2]i. However, one should 
be cautious with these values. First, the two datasets have a fre- 
quency difference of a factor of 10, which results in small formal 
error bars of the derived spectral index. One should be aware that 
the spectral index here could provide us with a trend in the spec- 
trum of HST-1 and the physical condition of the region, but not 
the absolute value. Second, we have found that the incomplete 
(m, v) coverage on short VLBA baselines has resulted in missing 
flux in the HST-1 region at A2 cm. Therefore, the spectral index 
value of HST-1 region is a lower limit, and the results suggest 
that HST-1 was optically thin during the flare in 2005. If assum- 
ing HST-1 was a flat-spectrum source, there would have been 
~90 mJy missing flux in this region in our 2 cm VLBA measure- 
ments. 

4.5. A blazar nature of HST-1 

The intensity of HST-1 reached a maximum in different wave- 
bands in 2005, and the light curves from the VLA A2 cm, VLBA 
/I20cm, NUV, and X-ray observations all show the same ten- 
dency ([Harris et al. 2009; Madrid 2009). Our Hght curve of HST- 
1 at VLBA /12 cm (Fig.|5]l is consistent with the other obser- 
vations. Although this trend could be taken as evidence that 
HST-1 is the source o f the HESS-detected TeV flare in 2005 
(lAharonian et al l l2006h . its VLBA radio properties that we have 
derived do not support it. These include (i) a very low com- 
pactness of the dominant emission, (ii) low brightness tempera- 
ture, (iii) sub-luminal motion, and (iv) possibly low optical depth 



across the feature at parsec-scales. Those are in contrast to the 
typical blazar core features, which tend to have flat or inverted 
spectral indices in cm-wave VLBA images. We observed that 
HST-1 has radio properties consistent with those commonly seen 
in jet components. TeV observatories cannot resolve the M87 
core and HST-1 separately. To probe the origin of TeV emission, 
coiTelating variability between different bands is a powerful 
tool. Recent examples of thi s approach were succ essfully shown 
for a sample o f blazars by Kovalev et al.l (12009 ). and for M 87 
bv lAcciari et al. (2009). Another approach is to apply physical 
models to multiwavelength observations to probe the TeV origin. 
There are models that suggest radio-TeV connections. For ex- 
ample, high energy flares could be generated from parsec-scale 
radio jets by inverse-Compton scattering of the photons and par- 
ticles emitting from the core, and Stawarz et al. (2006) used this 
approach to e xplain the TeV flare from M87 in 2005. However, 
lAcciari et al] t2009) reported about the VHE flare of M87 in 
2008, and suggested that there might be coiTelations with the ra- 
dio flare observed by VLBA at 43 GHz from the core. From our 
results, although we cannot fully discount, however, that HST- 
1 was the source of the TeV flare in 2005 based on our results 
alone, we do not favor the blazar nature for HST-1 suggested by 
iHarris et all ( 12008) . 

5. Summary 

With our VLBA A2cm data from 2000.06 to 2009.10, we have 
detected HST-1 during 2003.09-2007.10, which covered the 
multi-band flaring period of HST-1. The total flux density of 
HST-1 varied from 4-24 mJy in our detections; by comparing the 
images of VLBA A2 cm and A20 cm, we saw a steep spectrum 
with a > -0.8 in this region. The projected apparent speed of 
HST-1 derived from the brightness peak position is 0.61+0. 31c, 
which suggests a sub-luminal nature at VLBA A2 cm. 

Our results showed that HST-1 is extremely extended at 
parsec-scales, and has a steep spectrum. No compact feature 
with a brightness temperature higher than 9x10^ K is present 
in the A2 cm VLBA observations of this region of the M 87 jet, 
which implies that HST-1 does not have the properties of a stan- 
dard blazar core. Combining our findings, we do not find evi- 
dence of a blazar nature for HST-1 in the jet of M 87. 
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